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ABSTRACT 

An  algorithm  was  developed  to  schedule  courses  at  the  Naval  Guided 
Missiles  School  that  would  minimize  the  peak  demand  for  a  given  single 
resource. 

The  algorithm  will  provide  a  schedule  for  any  given  number  of  courses 
of  equal  length.  An  extension  to  the  algorithm  will  provide  a  schedule 
for  any  number  of  courses  of  unequal  length. 

The  algorithm  was  developed  in  two  phases.  The  first  phase  concerned 
finding  an  initial  schedule.  The  second  phase  concerned  reducing  the  max- 
imum demand  for  the  resource. 
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I.   INTRODUCTION 

A.   PROBLEM  MOTIVATION 

To  function  effectively,  a  weapons  system  as  vast  and  complex  as  the 
Fleet  Ballistic  Missile  system  (FBM) ,  requires  competent  and  highly  skilled 
technicians.   The  scheduling  of  the  required  training  courses  at  a  train- 
ing facility  is  a  vital  part  of  ensuring  the  FBM  system  is  manned  with  the 
necessary  expertise  to  operate  effectively. 

The  Office  of  the  Chief  of  Naval  Operations  (OPNAV)  determines  the  re- 
quirements for  technicians  to  man  the  FBM  force  based  upon  policy  objectives 
as  promulgated  by  the  Department  of  Defense.   The  Bureau  of  Naval  Personnel 
(BUPERS)  is  charged  with  the  responsibility  of  meeting  the  requirements 
specified  by  OPNAV.   BUPERS  devises  the  content  of  the  several  training 
curricula  necessary  to  fulfill  the  technician  requirements.   The  Naval 
Schools  Commands  are  assigned  the  responsibility  of  conducting  the  train- 
ing curricula.   Training  for  a  portion  of  the  FBM  personnel  is  conducted 
at  the  Naval  Guided  Missiles  School,  Dam  Neck,  Virginia.   At  Dam  Neck, 
courses  of  instruction  are  available  for  commanding  officers,  executive 
officers,  department  heads,  division  officers  and  the  following  technical 
ratings:   Fire  Control  Technician  (FT);  Electronics  Technician  (ET) ;  Mis- 
sile Technician  (MT) ;  and  Torpedoman  (TM) . 

The  training  cycle  for  FBM  enlisted  technicians  consists  of  recruit 
training,  "A"  school,  "B"  school,  submarine  school,  team  and  refresher 
training  and  operational  experience.   A  man  learns  the  basic  skills  of  his 
particular  rating  in  "A"  school.   The  "B"  school  provides  the  advanced 
technical  training  related  to  his  specific  job  in  the  FBM  system.   Team 


and  refresher  training  provide  periodic  review  of  basic  knowledge  and  in- 
corporate the  latest  technological  advices  of  the  FBM  system.  These  last 
three  areas  of  training  are  available  at  the  Naval  Guided  Missiles  School. 

The  FBM  system  is  presently  in  transition  from  utilization  of  the  POLARIS 
missile  to  the  POSEIDON  missile.   The  POSEIDON  missile  system  while  similar 
to  POLARIS  requires  technician  familiarization  training.   As  a  result,  the 
Guided  Missiles  School  is  faced  with  the  possibility  of  having  to  build  an 
additional  building  to  house  the  required  classrooms  and  laboratories  at  a 
projected  cost  of  $35.5  million. 

Given  the  resource  requirements,  such  as  number  of  instructors,  size  of 
classrooms,  size  of  the  laboratories,  for  the  convenings  in  the  training 
cycle,  a  schedule  to  minimize  the  peak  demands  for  these  resources  would 
keep  the  required  number  of  instructors,  classrooms  and  laboratories  as 
small  as  possible  and  thereby  keep  the  cost  of  the  building  as  low  as 
possible. 

Burgess  and  Killebrew  f_   1_/  present  an  heuristic  approach  to  scheduling 
project  activities  by  resource  leveling.   The  objective  they  use  is  the 
minimization  of  the  square  of  the  resource  requirement  in  each  time  seg- 
ment.  The  solution  seeks  to  keep  a  level  demand  for  the  resource  through- 
out the  scheduling  cycle.   If  the  conducting  of  a  course  by  the  training 
facility  is  defined  as  completing  a  project,  then  the  school  scheduling 
problem  may  be  viewed  as  a  multi-project  scheduling  problem. 

Willingham  j_   2_/  determined  the  optimum  number  of  convenings  for  a 
series  of  courses  in  a  one  year  time  period  utilizng  three  resources.   He 
then  presented  an  heuristic  method  to  schedule  the  convenings  into  the 
one  year  period.   He  attempted  to  balance  the  schedule  by  keeping  the  num- 
ber of  convenings  per  course  in  session  at  any  time  during  the  period 
constant. 


B.  DEFINITION  OF  THE  PROBLEM 

Given  the  required  student  input  into  the  FBM  system  as  determined  by 
higher  authority  and  the  required  teaching  resources  as  determined  by  the 
training  facility,  an  optimal  scheduling  process  is  necessary  to  make  eco- 
nomic use  of  the  resources. 

The  development  of  an  algorithm  that  will  minimize  the  peak  demands  of 
given  resources  is  required.   The  development  of  an  algorithm  that  will 
minimize  the  peak  demand  for  a  single  resource  is  the  first  step  towards 
such  an  algorithm. 

C.  SCOPE  OF  THE  STUDY 

The  problem  will  be  approached  by  developing  a  scheduling  algorithm  for 
a  single  resource  for  courses  of  equal  length  and  the  first  steps  of  a 
scheduling  algorithm  for  courses  of  unequal  length. 


II.   THE  SCHEDULING  ALGORITHM 

A.   PREVIEW 

The  central  idea  behind  the  algorithm  methodology  is  to  rearrange  the 
initial  schedule  in  order  to  reduce  the  maximum  demand  for  the  resource. 
The  algorithm  begins  by  scheduling  the  convenings  in  order  relative  to  their 
demand  for  the  resource;  those  with  the  greatest  demand  are  scheduled  first 
and  those  with  the  least  demand  last.   The  algorithm  then  attempts  to  inter- 
change convenings  scheduled  into  the  period  of  maximum  demand  and  conven- 
ings having  a  lesser  demand  for  the  resource  scheduled  into  any  other  time 
period,  thereby  reducing  the  maximum  demand. 

The  algorithm  presented  will  provide  a  schedule  for  courses  of  equal 
length  of  which  the  time  period  T  is  an  integer  multiple. 

Before  stating  the  algorithm  it  is  necessary  to  define  some  terms. 

Let  T  =  the  length  of  the  fixed  scheduling  period. 

I  =  number  of  different  courses  which  must  be  taught  during  T. 

N-£=  number  of  convenings  of  course  i,  (i=l  ,2  , .  .  .  .  ,1)  . 

tj[=  the  length  of  one  convening  of  course  i;  t^  ^T. 

c-£=  the  integer  resource  demand  coefficient  of  course  i. 
Clearly,  if  the  combined  lengths  of  all  required  convenings  does  not 
exceed  T,  that  is, 


/_  N.t. 


-CT, 
l 
then  the  courses  should  be  scheduled  sequentially  in  the  time  zero  to  T, 

/  0,T_/.   Thus  the  maximum  demand  for  the  resource  will  be  the  largest  of 

the  resource  demand  coefficients,  c. 


Assuming  the  problem  is  other  than  the  trivial  case  stated  above, 
a  scheduling  cycle  of  length  T  is  defined  as  a  fraction  (<T1.0)  of  a  sched- 
ule.  Hence  a  schedule  will  consist  of  at  least  two  cycles. 

Using  a  matrix  representation  for  the  problem,  an  mxn  matrix  lo.  .1  may 


H 


be  defined  where  o..- .  denotes  the  resource  demand  coefficient,  c,  .  of  the 
convening  scheduled  into  the  j    segment  of  the  i   scheduling  cycle. 

The  dimension  of  the  matrix,  mxn,  is  determined  by  the  algorithm;  n  is 
equal  to  the  number  of  time  segments  in  the  time  period  T,  and  ra  is  twice 
the  number  cf  scheduling  cycles  needed  to  complete  the  schedule  minus  one. 

The  problem  to  be  solved  by  the  algorithm  may  be  stated  as,  find  that 
schedule  of  course  convenings  which  minimizes 


=  max  2~o-.    , 


J    i 
the  largest  resource  demand  to  occur  at  any  time  during  the  period  T. 

B.   THE  ALGORITHM  FOR  COURSES  OF  EQUAL  LENGTH 

The  algorithm  begins  by  determining  an  initial  schedule.   An  attempt 
is  then  made  to  decrease  the  maximum  demand  for  the  resource  by  reschedul- 
ing any  combination  of  courses.   The  reschedule  is  first  attempted  between 
the  period  of  maximum  demand  and  the  period  of  minimum  demand.   If  this  is 
not  possible,  the  algorithm  then  tries  to  reschedule  courses  between  the 
period  of  maximum  demand  and  any  other  period. 

After  the  initial  schedule  is  obtained,  the  algorithm  compares  each 
single  course  scheduled  into  the  period  of  maximum  demand  with  each  single 
course  scheduled  into  the  period  of  minimum  demand  in  an  attempt  to  find 
two  courses  whose  interchange  will  reduce  the  maximum  demand  for  the  re- 
source.  If  no  single  course  can  be  found  in  the  minimum  period  to  be  used 
in  an  interchange,  the  algorithm  successively  compares  each  single  course 


in  the  maximum  period  with  all  pairs,  triplets,  etc.,  of  courses  in  the 
minimum  period.   If  the  total  demand  for  these  combinations  of  courses  will 
not  reduce  the  maximum  demand  when  interchanged  with  a  single  course  of  the 
maximum  period,  the  algorithm  compares  first  pairs,  then  triplets  and  so  on 
of  courses  in  the  maximum  period  with  combinations  of  courses  in  the  mini- 
mum period.   If  after  all  possible  combinations  of  courses  in  the  maximum 
period  are  compared  with  all  combinations  of  courses  in  the  minimum  period 
and  no  interchange  of  courses  can  be  made,  the  algorithm  selects  another  time 
period  in  place  of  the  minimum  period  and  starts  the  process  again.   After 
all  time  periods  have  been  investigated  and  fail  to  produce  an  interchange, 
the  algorithm  terminates. 

1 .   Statement  of  the  Algorithm 

The  number  of  time  segments  in  the  time  period  T  can  be  calculated 


by 


TS  =   T 


where 


TS  =  the  number  of  time  segments. 
t  =  the  length  of  one  convening  of  any  course. 


1.   There  is  a  resource  demand  coefficient  c^  for  each  convening  of 
course  i.   Hence  there  are  a  total  of  2_  N •  coefficients.   Arrange  these 
demand  coefficients  into  a  non- increasing  sequence  \sA    where  j  =  1 ,2 , .  .  .  ,£_N- , 
Compute  P  where 


-Ei 


TS 


If  P  -  H^    i 


s  an  integer,  then  add  (P-l)TS  zero  elements  to  the  end 


TS 
of  the  sequence.   If  P  is  not  an  integer,  then  add  (K-P)TS  +  (K-1)TS=(2K-P-1)TS 


zero  elements  to  the  end  of  the  sequence  where  K  is  the  smallest  integer 
larger  than  P.   Schedule  the  first  cycle  of  convenings  by  placing 

Si ,S2> sTq  into  the  (m,l) , (m, 2) ,...., (m,TS)  positions;  e.g.,  the  bottom 

row  of  the  matrix.   The  second  cycle  is  obtained  by  scheduling  backwards 
over  the  period  jO,TJ  and  adding  the  (m-1)    row  to  the  matrix.   Thus,  this 
row  consists  of  s  Tq  ,s„T„_,  ,...., s„,  .  The  schedule  and  the  matrix  are 
completed  by  continuing  to  schedule  forward  on  the  odd  numbered  cycles  and 
backwards  on  the  even  numbered  cycles  until  the  sequence  \s-l   is  exhausted 

2.  Compute  the  value  of  o  where 

0  =  max  2__  o-y  -  min  / O-jj  . 

J    i         J    i 

If  0  =  0  or  1 ,  then  stop;  the  solution  is  optimal.   Otherwise  go  on  to 
step  3. 

3.  Designate  the  columns  of  maximum  and  minimum  sums  from  the  matrix  as 

o  and  o   ,  respectively.   In  case  of  ties,  select  a  column  at  random.   Let 
ou  and  ov  represent  some  single  element  of  ou  and  ov   .   Set  x  =  o   and 
go  to  step  4. 


4.   Set  y  =  max  \  o-^v  <£x  (   .   Interchange  the  elements  of  o   and  o   cor- 


\°iv<xl 
responding  to  x  and  y  if  x  ~(j<^y  <^  x  and  go  to  step  10.   Otherwise,  reduce 

x  to  the  next  largest  element  value  in  o   and 

a   if  x^>0,  repeat  step  4. 

b.   if  x  =  0,  set  x  =  o   ,y,  =  max   \  o^v /x  and  go  to  step  5. 

5.   Denote  the  element  corresponding  to  y,  as  o   . 

Set  y2  =  max  io±v  j   i<p,  0<oiv^c  -  yA      . 

a.   If  an  element  of  o   can  be  found  to  use  as  y?  then  go  to  step  6. 


10 


b.   If  no  element  of  o   can  be  found  to  use  as  y   then  reduce  y,  to  the 
value  of  the  element  in  o   which  is  the  next  smaller  than  its  current  value 
and  again  try  to  find  a  yo  value.   If  no  y0  value  is  found  even  when  y,  is 
reduced  to  the  second  smallest  non-zero  element  in  o  ,  go  to  step  9. 

6.  Interchange  the  elements  of  o"v  corresponding  to  yj  and  y2  with  the 
element  of  ~ou   corresponding  to  x  and  any  zero  element  of  o"  if 

x  -  ^  <^yi  +  yo  -^x  ar-^  g°  to  step  10.   Otherwise,  go  to  step  7. 

7.  Denote  the  element  corresponding  to  y2  as  o    . 

Set     y3  =  max  f  o.v    i  <  q  ,  0  <\>±v  <x  -  yj_    -  yJ  < 

a.  If  an  element  of  o"  can  be  found  to  use  as  y0  then  set  k  =  3  and 

v  J  3 

go  to  step  8. 

b.  Otherwise,  reduce  y2  to  the  value  of  the  next  smaller  positive  ov 
element  if  one  exists  and  return  to  step  6.   If  no  positive  element  exists 
to  use  as  y  return  to  step  5b. 

8.  Continue  the  process  begun  in  steps  6  and  7  of  trying  to  interchange 
groups  of  elements  corresponding  to  yi  ,  yo,....,yk  from  o  with  the  element 
corresponding  to  x  and  k-1  zero  elements  in  ou.   For  example,  search  next 
for  a  y4  satisfying  y4  =  max  j  oiv   i <r ,o  <o±v  <x-yi~y2-y3 j  .   If  no  y 4 
exists  then  reduce  y~  and  again  attempt  an  interchange  with  y, ,  y- ,  and  y„ 
if  y3>0.  If  y3  =0,  go  to  step  7b. 

9.  Decrease  x  to  the  value  of  the  next  largest  element  of  o"  .  If  x^O, 
set  y^  =  max  )  o^v  -^ x  j  and  return  to  step  5.  Otherwise,  set  x  =  o  + 
o      and  go  to  step  11. 

10.  Rearrange  the  elements  of  o"  and  o"  into  a  non-increasing  sequence 
beginning  with  the  bottom  element  of  each  column.   Go  to  step  2. 


11 


11.  Set  y  -  max  j  o  •  <^x  j  .   Interchange  the  elements  of  o"  and  o"  cor- 
responding to  the  elements  used  in  the  x  and  y  values  if  x  -  J  <^y^.x  and 
go  on  to  step  10.   Otherwise,  reduce  x  as  follows: 

current  x  value  =  opu  +  op_2  u  ; 
new  x  value  =  o^  ^   +  op_2)U. 

a.  If  x>0,  repeat  step  11. 

b.  if  x  =  0,  set  x  =  o   +  o   i_  n^d   proceed  as  in  steps  5  through  9 

in  an  attempt  to  exchange  a  pair  of  elements  of  o"  with  two  or  more  elements 
from  Oy.   If  no  interchanges  are  possible  go  to  step  12. 

12.  Successively  increase  the  number  of  elements  in  o"  used  in  the  computa- 
tion of  x  by  one  more  element  and  continue  the  process  of  first  comparing 
these  elements  with  single  elements  of  o  and  then  with  pairs,  triplets,  etc. 

If  no  interchange  between  o"  and  ~o     are  possible  then  define  o"  as  the  col- 
6  u      v     ^  v 

umn  having  the  next  smallest  sum  and  return  to  step  2. 

2.   Example  1 

This  example  will  demonstrate  the  forward  and  backwards  scheduling 
method  for  obtaining  an  initial  schedule. 

Let  us  suppose  a  training  facility  must  offer  five  courses  during  a 
ten-month  training  cycle  according  to  the  following  table: 


12 


Table  1.   Problem  Description  for  the  Example  of  Case  I 


Training 
Course  i 

1 

2 

3 

4 

5 


N 


l 
6 

9 

1 

8 

10 


2 

9 

8 
3 

4 


where 


N-  =  the  required  number  of  convenings  of  course  i. 

tj[  =  the  length  in  months  of  one  convening  of  course  i. 

c^  =  the  required  laboratory  space,  in  square  feet,  for 
one  convening  of  course  i. 


The  resource  demand  coefficients  c-  are  arranged  into  a  non-increasing  se- 
quence: ^9,9,9,9,9,9,9,9,9,8,4,4,4,4,4,4,4,4,4,4,3,3,3,3,3,3,3,3,2,2,2,2,2,2, 

0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0!  . 

The  total  time  period,  T,  is  ten  months.   TS  is  determined  to  be  10.   There- 
fore the  segment  length  must  be  one  month.   The  matrix  will  be  a  mxlO  matrix. 
The  first  scheduling  cycle  is  completed  by  scheduling  the  first  ten  ele- 
ments of  the  sequence.    / 

999999999 

The  second  cycle  is  completed  by  scheduling  the  next  ten  elements  of  the 
sequence  backwards  over  the  period. 


4444444444 
999999999 
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The   third   cycle   is   completed   by  forward   scheduling   elements   twenty-one 
through   thirty.  /  a 


3333333322 

4444444444 


999999999 


The   fourth  cycle   is   completed  by  scheduling   the  next    ten  elements 


0000002222 
3333333322 
4444444444 

9999999998 


?••<- 


16  16  16  16  16  16  18  18  17  16 


The  matrix  is  completed  by  adding  the  zero  elements. 


Column  seven  is  designated  o  .   Column  six  is  designated  o  . 

to        u  b  V) 


0000000000 
0000000000 
0000000000 
0000002222 
3333333322 
4444444444 
'9  999999998 

Ho..      =   16  16  16  16  16  16  18  18  17  16 
i  1J 

0  =18-16 

=  2.   The  value  for  x  =  o   is  9  as  required  by  step  4  of  the  algorithm. 

For  x  =  9,  a  value  for  y  satisfying  x  -  (j  -^yCx  cannot  be  found. 

The  value  of  x  is  reduced  to  4,  the  next  largest  element  in  o"   ,  as 

specified  by  step  4  and  the  search  for  y  is  repeated.   A  value  of  y  =  3 

satisfies  the  criterion.   Interchanging  the  elements  corresponding  to  x  and 

y,  rearranging  the  columns  required  by  step  10  not  being  necessary,  results 

in  the  schedule: 
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0000000000 
0000000000 
0000000000 
0000002222 
3333343322 
4444443444 
999999999  8, 


£°ij  = 


r 


16  16  16  16  16  17  17  18  17  16 
Returning  to  step  2,  column  eight  is  ou,  column  ten  is  o"  and  d    =  1< 
16  =  2.   Proceeding  as  before  in  the  search  for  an  acceptable  value  for  y, 
step  4  is  repeated  until  x  has  been  reduced  to  3.   A  value  for  y  =  2  quali- 
fies to  be  interchanged  with  x.   Rearranging  the  matrix  as  required  by  step 
10  produces  the  schedule: 


l 


0  0  0  0  0 
00000000 
00000000 
00000022 
3  3  3  3  3 

4  4 

9  9 


0  0  0  0  0  \ 


0  0 
0  0 


4 
\  9 


4  4 
9  9 


=  16  16  16  16  16  17  17  17  17  17 


The  schedule  is  optimal  beca 


use()= 


17  -  16 


1. 


3.   Example  2 

Let  us  suppose  the  initial  schedule  for  a  particular  problem  is 
given  below. 


27  25 


Column  one  is  designated  o  .   Column  two  is  designated  o  . 


U- 


25  -  2, 


The  value  for  x  =  0   is  13  as  required  by  step  4  of  the  algorithm.   For 


mu 
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x  =  13,  a  value  for  y  of  step  4  such  that  x  -  (\  ^ly  <Cx  cannot  be  found. 

The  value  for  x  is  reduced  to  7,  the  next  largest  element  in  "o  .   Again 
no  value  for  y  will  satisfy  the  criterion  x  -  fl  <Cy  <Cx.   This  is  also  true  for 
x  =  5  and  x  =  2.   The  value  for  x  is  reduced  to  zero  and  the  algorithm  looks 
for  a  pair  of  elements  in  o"v   to  interchange  with  x. 

The  value  for  x  is  again  set  to  o    =  13.   The  largest  element  of  b~ 
less  than  x,  designated  y,,  is  10.   The  second  element  of  the  pair,  y2,  must 
satisfy  10  +  y2\13  with  y  >0.   No  such  y   can  be  found.   Therefore,  the 
value  of  y,  is  reduced  to  9  as  specified  by  step  5  and  the  search  for  y~  is 
repeated.   This  time  a  value  of  3  is  found  to  use  as  y2-   The  pair  (9,3)  from 
the  second  column  can  therefore  be  interchanged  with  the  element  13  and  any 
of  the  0  elements  in  the  first  column.   Rearrangement  of  the  elements  ac- 
cording to  step  10  results  in  the  following: 


Lo±i 


0 

o\ 

'   0 

0  \ 

2 

0 

3 

0    ] 

i  5 

3    ) 

\7 

10  j 

13  / 

*9 

26  26 


6  =  0  u 


i 


Since  0=0  upon  returning  to  step  2,  the  solution  is  optimal. 

C.   THE  ALGORITHM  EXTENDED  TO  COURSES  OF  UNEQUAL  LENGTH 

When  the  scheduling  problem  is  extended  to  courses  of  unequal  length, 
it  becomes  more  difficult.  Each  convening  will  now  span  one  or  more  time 
segments  depending  upon  the  individual  course  length.   It  can  no  longer 
be  guaranteed  that  a  scheduling  cycle  will  be  complete.   There  may  be  gaps, 
due  to  the  unequal  course  lengths,  where  no  convening  may  be  scheduled. 
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The  solution  method  becomes  more  complicated  in  that  interchanging 
of  convenings  to  reduce  the  maximum  demand  for  the  resource  must  make  allow- 
ance for  the  unequal  course  lengths.   The  algorithm  must  be  modified  such 
that  reducing  the  maximum  demand  will  consist  of  a  series  of  one  or  more 
convenings  being  interchanged  with  another  series  of  convenings  spanning 
the  same  number  of  time  segments. 

The  algorithm  presented  in  this  section  considers  only  interchanging 
individual  convenings.   Clearly,  interchanging  series  of  convenings  instead 
of  individual  convenings  becomes  especially  difficult  when  the  algorithm 
must  consider  pairs,  triplets,  etc.,  of  convenings  to  interchange  and  will 
be  left  as  the  subject  of  a  future  study. 

1 .   Statement  of  the  Algorithm 

The  number  of  time  segments  in  the  time  period  T  can  be  calculated 

by  redefining  the  variables  TS  and  t. 

Let  TS  =  T 

t 

where 

TS  =  the  number  of  time  segments 

t  =  the  greatest  common  divisor  for  all  course 
lengths  and  T. 

Due  to  the  unequal  lengths  of  the  courses,  some  convenings  may  span 

more  than  one  time  segment.   The  number  of  time  segments  spanned  by  each 

convening  is  the  ratio 

TS±  =  ti 
t 

where 

TS.  =  number  of  time  segments  spanned  by  one  convening 
of  course  i 

t^  =  length  of  one  convening  of  course  i. 
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1.  There  is  a  demand  coefficient  for  each  convening  of  course  i.   Each  con- 
vening may  be  represented  by  a  group  of  identical  coefficients  indicating  the 
number  of  time  segments  spanned  by  the  convening;  e.g.,  (c.,c.,c.)   if  course 
i  spans  three  segments.   There  are  {_   N-  coefficient  groups.   A  convening  of 
only  one  time  segment  in  length  will  be  represented  by  a  single  integer  with- 
out parentheses.   Arrange  these  demand  coefficients  into  a  non-increasing 

sequence  of  groups  \   (s  . ,s  ......  ,s;  j  )(  where  j  =  1,2,...,  LN.  .   Begin 

12        TS  • 
scheduling  the  first  cycle  of  convenings  by  placing  s1  ,s.  ,...., s1  1)  into 

the  (m,l),  (m, 2)  ,....,  (m,TS,  )  positions  of  the  matrix.    Com- 
plete the  first  cycle  by  placing  groups  into  the  (m,TS^  +  1),  (m,TSj  +2),..., 
(m,TS)  position  until  the  (m,TS)  position  is  filled  or  a  complete  group  can- 
not be  scheduled  into  the  time  remaining  due  to  its  length.   In  this  latter 
situation  the  bottom  matrix  row  is  completed  with  zero-element  groups.   The 
(m-1)  row  is  then  scheduled  using  the  backwards  pass  and  adding  zeros  when 
needed  to  fill  up  the  row.   The  schedule  and  the  matrix  are  completed  by 

continuing  to  schedule  forward  on  the  odd-numbered  cycles  and  backwards  on 

!1    2        TS •   / 
(s   ,s   , . . . . ,s  J  )  is 
j    j        j 
J 
exhausted. 

2.  Compute  the  value  of  0  where 

0    =   max  / o  .  .  -  min  /    o  .  . . 

J    1        J    1 

If   0  =  0  or  1 ,  then  stop;  the  solution  is  optimal.   Otherwise  go  to  step  3. 


3.   Designate  the  column  of  maximum  sum  from  the  matrix 

as  o  .   In  case  of  a  tie,  select  a  column  at  random.   Let  o   represent  some 

u  '  »  u   ^ 

single  element  of  o"  .   Define  the  group  containing  o  as  (o   ,0^  , .  .  .  .  ,o-1'^w) 

U  or  °u  uuu 

where  w  is  the  position  number  of  the  course  for  which  o  is  the  resource 

u 


demand  coefficient  in  the  demand  coefficient  sequence  of  step  1.   Let  o 


v 
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represent  any  other  single  element  of  the  matrix  not  contained  in  ou.   De- 
fine the  group  containing  ov  as  (oj  ,aj  ,...,oJSy  )  where  y  is  the  position 
number  of  the  course  for  which  ov  is  the  resource  demand  coefficient  in  the 
demand  coefficient  sequence.   Let   \(s^    ,s^   >--->sk   >J  ijj  be  a  £ubse_ 
quence  of  P  demand  coefficient  groups  starting  with  the  (i,j)  position  of 

the  matrix  and  spanning  \    <TS.time  segments  with  \    -  2_  TSk-   Each  sub- 

k=l 
sequence  will  consist  of  a  portion  (  <1)  of  row  i.   Set  ou  =  the  maximum 

valued  non-zero  integer  element  in  ou- 

Search  the  matrix  for  an  element  ov  which  satisfies  the  following  con- 


ditions 


i.   o   -c^<ott<o   where  o<=  )_  °iu  ~  /Z° 
U         V     U  .    xu-         ^-7- 


IV 

V 


ii.    (oj  ,o|  ,...,oTSy  )  £   |^(sk  ,Sk  ,...,SkTSk)]  ±n 

for  some  i,j,  )  having  \    =  TSw 

*"•    <«J  .•* Efr   >  fr  [(si  ,sk skSk^« 

has  sr  ^0u  for  all  r.                         12 
If  such  an  element  can  be  found,  interchange  the  group  (ou  ,ou  , , 

o^w  )  with  the  subsequence   j(s£  ,s£  ,...,sJ;Sk  )|  . j j  in  their  respec- 
tive columns  and  return  to  step  2.   If  no  interchange  is  possible  set  ou 
equal  to  the  next  largest  integer  in  o~u  and  repeat  the  search.   If  no  in- 
terchange is  possible  for  ou  being  any  non-zero  element  of  o"u,  then  terminate 
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2.   Example  1 

Let  us  suppose  the  training  facility  must  now  offer  courses  of  vari- 
ous lengths  as  given  in  Table  II.   Again  it  is  desired  to  minimize  the  peak 
demand  for  the  resource. 

Table  II:   Problem  Description  for  the  Example  of  Case  II 


Ci 


Training 
Course    ± 

»i 

H 

1 

4 

l 

2 

6 

3 

3 

4 

4 

4 

8 

2 

5 

4 

1 

1 

2 
3 
4 
6 


where 


c  -  the  number  of  classrooms  required  for  the  student  enrollement 


i 

of  course  l. 


The  training  cycle,  or  time  period  T,  is  specified  as  one  year.   The  great- 
est common  divisor  of  all  course  lengths  and  the  time  period  is  one  month. 
Therefore  TS  is  equal  to  12  time  segments.   The  values  for  TSi  can  be  cal- 
culated as  given  in  Table  III. 

Table  III:   Time  Segments  per  Course  for  the  Example  of  Case  II 

Training 

Course  i  TS-j_ 

1  1 

2  3 

e 

3  4 

4  2 

5  1 
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The  sequence  of  resource  demand  coefficients  is:  j  6,6,6,6,(4  4)  ,  (4  4) 
(4  4), (4  4),  (4  4),  (4  4),  (4  4),  (4  4),  (3  3  3  3),  (3  3  3  3) ,  (3  3  3  3), 
(3  3  3  3),  (2  2  2),  (2  2  2),  (2  2  2),  (2  2  2),  (2  2  2),  (2  2  2),  1,1,1,1  ^ 

There  are  twelve  time  segments  in  the  schedule,  therefore  the  segment 
length  is  one  month.   The  matrix  will  be  an  mxl2  matrix. 

The  initial  solution  is  obtained  using  the  forward  and  backward  method 
discussed  in  Case  I. 

'(2  2  2)  (2  2  2)  1  1   1   1   0  0 

(2  2  2)  (2  2  2)  (2  2   2)  (2   2  2) 

(3  3  3  3)  (3  3  3  3)  (3  3  3  3) 

\(3  3  3  3)  (4  4)(4  4)  (4  4)(4  4) 

\6  6  6  6  (4  4)(4  4)(4  4)(4  4), 


£  o 


13 


16  16  16  16  .  15  15  14.  14.  14  14  13  13 


From  step  2  of  the  algorithm 

I 


we  see 


r     y~ 

0    =  max  i—r-  o  •  •  - 


i   iJ 


mm  T"  o-h  =  16-13  =  3.   Therefore  some  rearrangement  of  the  matrix  may 
be  necessary.   Let  column  one  be  designated  as  the  maximum  column  and  set 
o  =  6,  the  largest  element  of  column  one.   For  o  =  6,  condition  i  of  step 
4  of  the  algorithm  is  violated  for  all  possible  values  of  o  . 

For  ou  =  3,  ov  =  2  will  satisfy  condition  i.   However,  condition  ii  is 
violated.   Reducing  o   to  2  and  setting  o   =  1  produce  o(  =  2  and  the  con- 
ditions are  satisfied.   The  group  (2  2  2)  and  sub-sequence   j 1,1,1 > 


are 


interchanged. 


1J 


(2  2  2)  (2   2  2)  (2  2   2)  10  0 

1  1  1  (2   2   2)  (2  2   2)  (2   2  2) 

(3  3  3  3)(3  3  3  3)  (3  3  3  3) 

(3  3  3  3)  (4   4)  (4  4)  (4  4)  (4  4) 

6  6  6  6  (4  4)  (4  4)(4  4)(4  4) 

15  15  15  16  15  15  15  15  15  14  13  13 


Returning  to  step  2,  q  is  again  computed  to  be  16-13  =  3.  Column  four 
is  now  the  maximum  column.  Set  o  =  6,  the  largest  element  of  column  four, 
and  search  the  matrix  for  values  of  ov.   As  before,  an  acceptable  value  for 
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o   cannot  be  found  until  ou  has  been  reduced  to  2.   Setting  o   =  1,  the 
group  (2  2  2)  is  interchanged  with  the  sub-sequence  £l,0,0j 

(2   2   2)  (2   2  2)  (2  2   2)  (2   2  2)' 

111  1   0  0  (2  2   2)  (2   2  2) 

(3  3  3  3)  (3  3  3  3)  (3  3  3  3) 

(3   3  3  3) (4  4)  (4  4) (4  4)(4  4) 

,666  6  (4  4)  (4  4)  (4  4)  (4  4)/ 

/    o . .   =   15  15  15   15  13  13  15  15  15   15  15  15 

X 

Trying  to  reduce  the  maximum  value  of  15  in  column  one  and  two  results 
in  two  interchanges  when  o  equals  1. 

(2   2   2)  (2   2   2)  (2   2   2)  (2   2   2) 

001  111  222    222 

(3  3  3  3)  (3  3  3  3)  (3   3   3  3) 

(3  3   3  3)  (4  4)  (4  4)  (4   4)  (4  4) 

6  6  6  6  (4  4)  (4  4)  (4   4)  (4  4) 

/_   o..  =  14  14  15  15  14  14  15  15  15  15  15  15 

1 

The  solution  is  optimal  because  (j  =  15-14  =  1. 

It  is  interesting  to  note  that  the  last  two  interchanges  are  not  neces- 
sary to  ensure  optimality.   Because  the  number  of  columns  (10)  having  the 
maximum  sum  of  15  is  greater  than  the  number  of  columns  (2)  having  the  mini- 
mum sum  of  13  and  it  would  never  be  possible  to  reduce  the  max  L„    o^  • 

i 
below  15. 

D.   COMMENTS  ON  OPTIMALITY 

1 .   The  Algorithm 

There  are  three  conditions  under  which  the  algorithm  will  terminate, 
The  first  two,  q    =  0  or  q    =  1 ,  will  produce  an  optimal  schedule.   There 
is  no  guarantee  of  optimality  if  the  algorithm  terminates  with  q   ^>  1 . 

If  O  =0,  then  the  maximum  column  sum  and  the  minimum  column  sum  are 

equal.   Clearly  all  column  sums  are  equal  and  the  z  value  is  constant  over 

the  time  period  T.   Therefore  the  interchanging  of  elements  between  time 

segments  cannot  reduce  this  value  and  the  solution  is  optimal. 
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If  O— i>  then  the  maximum  column  sum  and  the  minimum  column  sum  differ 
by  only  one  unit.   Since  the  elements  of  the  matrix  are  non-negative  integers, 
interchanging  elements  between  the  maximum  and  minimum  columns  will  have  one 
of  three  effects.   The  value  of  z  will  not  change  if  the  elements  are  equal. 
If  the  elements  differ  by  one  unit  and  o   <C°U»  tne  maximum  and  minimum  col- 
umns will  become  the  minimum  and  maximum  columns,  respectively.   If  they  dif- 
fer by  one  unit  and  o  >  o   then  the  maximum  value  would  increase.   The 
maximum  value  will  also  increase  if  they  differ  by  more  than  one  unit.   In 
no  case  will  z  be  reduced  and  the  solution  is  optimal. 

For  A^l,  the  maximum  column  and  the  column  sums  of  all  other  columns 
differ  by  an  amount  less  than  the  combined  difference  between  any  combination 
of  elements  in  o"u  and  any  other  column.   It  is  possible  that  some  intermediate 
interchange  between  two  o   columns  may  allow  a  subsequent  interchange  to  take 
place  between  some  o"v  and  the  current  o"  column.   Therefore  no  general  guar- 
antee of  optimality  can  be  provided. 

2.   The  Algorithm  Extended 

The  algorithm  terminates  under  one  of  three  conditions.   If  0   = 
0  or  1 ,  the  previous  discussion  on  optimality  is  germane. 

If  $>1,  the  algorithm  attempts  to  reschedule  any  group  of  elements  in 
the  maximum  column  that  will  reduce  the  z  value.   There  is  no  provision  in 
the  algorithm  to  reschedule  combinations  of  groups.   For  this  reason  there 
can  be  no  guarantee  of  optimality. 
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III.   DISCUSSION  AND  EXTENSIONS 

A.  DISCUSSION  OF  THE  ALGORITHM 

To  obtain  an  initial  solution  to  the  algorithm  and  the  extension,  the 
sequence  of  convenings  is  placed  into  the  matrix  by  sched- 
uling forward  on  the  odd-numbered  cycles  and  backwards  on  the  even-numbered 
cycles.   While  this  is  not  the  only  way  to  initially  complete  the  matrix,  it 
is  designated  to  accommodate  the  user.   Because  the  algorithm  is  not  de- 
pendent upon  the  method  of  obtaining  an  initial  solution,  the  forward  and 
backwards  scheduling  seems  to  be  an  easy  way  to  get  started. 

B.  ADDITIONAL  OPTIMALITY  CONDITIONS 

There  appear  to  be  tests  other  than  Q    =  0  or  1  for  determining,  at  any 
time  during  the  solution  process,  if  optimality  has  been  obtained.   One  such 
test  involves  problems  having  multiple  columns  of  maximum  sum  and  multiple 
columns  of  minimum  sum.   For  example,  in  the  case  of  o    =  2,  if  the  number 
of  multiple  columns  of  maximum  sum  is  greater  than  the  number  of  multiple 
columns  of  minimum  sum,  the  maximum  demand  cannot  be  reduced  and  the  solu- 
tion is  immediately  optimal.   For  (\  >3,  the  relationship  of  the  number  of 
multiple  columns  of  maximum  sum  to  the  number  of  multiple  columns  of  mini- 
mum sum,  for  which  the  solution  is  immediately  optimal,  is  not  easily 
determined. 

Further  study  of  optimality  conditions  should  also  involve  interchanges 
between  columns  other  than  the  maximum  and  minimum  columns.   Such  a  study 
may  require  an  analysis  of  a  large  number  of  randomly  generated  examples. 
An  analysis  of  this  type  would  also  be  beneficial  in  appraising  how  well 
the  algorithm  obtains  optimal  solutions. 
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C.   EXTENSION  TO  MORE  THAN  ONE  RESOURCE 

The  scheduling  problem  where  more  than  one  resource  is  involved  in  com- 
plicated by  the  fact  that  minimization  of  each  ignoring  the  others  will  un- 
doubtedly result  in  completely  different  schedules,  any  of  which  does  not 
minimize  all  resources.   One  approach  would  be  to  create  a  combined  resource 
by  weighting  the  individual  resources  according  to  their  importance  and 
adding  them  together.   The  algorithm  described  in  this  thesis  could  then  be 
used  to  find  a  good  schedule. 

The  more  common  problem  is  to  minimize  the  demands  for  one  unrestricted 
resource  subject  to  restrictions  on  the  maximum  amount  to  resources  availa- 
ble for  the  remainder  of  the  resources.   The  preceding  algorithms  would  have 
to  be  modified  to  generate  an  initial  feasible  schedule  and  then  to  insure 
that  any  column  interchanges  which  take  place  do  not  violate  feasibility. 
The  development  of  such  modifications  should  begin  using  the  first  model 
(all  courses  of  equal  length)  since  it  provides  the  simplest  structure  for 
checking  feasibility  of  proposed  interchanges. 
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